Mott Insulating Ground State and its Proximity to Spin-Orbit Insulators in Na 2 Ir0 3 
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We present an anti-ferromagnetically ordered ground state of Na2lr03 based on density-functional-theory 
calculations including both spin-orbit coupling and on-site Coulomb interaction U. We show that the splitting of 
e' g doublet states by the strong spin-orbit coupling is mainly responsible for the intriguing nature of its insulating 
gap and magnetic ground state. Due to its proximity to the spin-orbit insulator phase, the magnetic ordering as 
obtained with finite U is found to exhibit a strong in-plane anisotropy. The phase diagram of Na2lr03 suggests 
a possible interplay between spin-orbit insulator and Mott anti-ferromagnetic insulator phases. 



PACS numbers: 71.70.Ej, 75.30.Kz, 71.20.-b, 75.30.Gw 

Recently, the role of spin-orbit coupling (SOC) has at- 
tracted great attention in many fields of condensed mat- 
ter physics. In multiferroic materials, for example, SOC 
combined with a large electron-lattice interaction has been 
suggested to be responsible for the multiferroic behavior 
which exhibit both non-collinear magnetic ordering and lat- 
tice polarization^. SOC is also indispensable to anomalous 
Hall and spin Hall effects where Hall and spin Hall currents 
are generated by an external electric field, respectively^™ In 
particular, the quantum spin Hall effect has led to the notion of 
topological insulators, new states of quantum matter 6 7 . While 
they have bulk energy gaps generated by the SOC, topological 
insulators are characterized by the presence of gapless surface 
states which are protected by time-reversal symmetry 8 . 

Another manifestation of strong SOC combined with on- 
site Coulomb interactions is the j e ff=l/2 Mott insulator dis- 
covered in Sr2lr04, one of the 5d transition-metal oxide d^HD. 
The novel spin-orbit integrated state with j eff =l/2 arises from 
the combined action of both strong SOC and intermediate on- 
site Coulomb interactions within the Ir 5d ti g manifold. In 
addition, there has been a theoretical proposal on the room 
temperature quantum spin Hall effect in Na2ir03 based on the 
j e ff=l/2 physics 11 , where the honeycomb lattice consisting of 
edge- shared Ir06 octahedra in each Ir-0 layer was considered 
to be an ideal realization of the Kane-Mele model, where hop- 
ping integrals between the j eff =l/2 states at the Fermi level 
was assumed to be an essential ingredient for the quantum 
spin Hall effect 8 12 . Since the crystal structure and local en- 
vironment of Ir atoms in Na2ir03 are different from those of 
Sr2lr04, however, it is necessary to clarify the electronic and 
magnetic structures of the Ir 5d manifold in this Na2ir03 com- 
pound with hexagonal lattice. 

In this paper, we present novel electronic structure and 
magnetic properties of Na2lr03 by carrying out density- 
functional-theory (DFT) calculations including both SOC 
Aso and on-site Coulomb interaction. We observe that a new 
form of spin-orbit coupled states emerges from the e' g dou- 
blet states near the Fermi level (Ep) and determines the in- 
triguing nature of its insulating gap. With an effective on-site 
Coulomb interaction parameter U = 2.0 eV, the ground state 
of Na2lr03 is found to be an antiferromagnetic (AFM) insu- 
lator with the ordered moments lying down within the honey- 
comb lattice of Ir atoms. The large splitting of the e' g doublet 



by the strong SOC is related to the strong in-plane anisotropy 
of magnetic ordering. Considering the role of SOC, we pro- 
pose a phase diagram in the Aso~^ parameter space which 
features a phase boundary between AFM Mott insulators and 
SO insulators. By estimating the exchange couplings between 
neighboring Ir atoms, we suggest a possible frustration of 
magnetic ordering in its ground state, which is consistent with 
a recent experiments 

In order to examine the effects of both SOC and on-site 
Coulomb interaction on the electronic structure of Na2ir03, 
it is necessary to treat both SOC and U on an equal foot- 
ing in the description of Ir 5d states. To identify the role of 
each contribution as well as the interplay between them, we 
carried out DFT calculations within the local-density approxi- 
mation (LDA), LDA including SOC (LDA+SO), and LDA+/7 
including SOC (LDA+/7 +SO) respectively. For the calcula- 
tions, we used the DFT code, OpenMxP, based on the linear- 
combination-of-pseudo-atomic-orbitals method 15 , where both 
the LDA+L 7 ^ method 16 and the SOC contribution were in- 
cluded via a relativistic j -dependent pseudo-potential scheme 
in the non-collinear DFT formalism. Double valence and sin- 
gle polarization orbitals were used as basis sets, which were 
generated by a confinement potential scheme with cutoff radii 
of 7.0, 7.0 and 5.0 a.u. for Na, Ir, and O atoms respectively. 
We used a (14 x 14 x 14) k-point grid for the k-space integra- 
tion. 

Up to our knowledge there is no crystal structure data for 
Na2ir03 published yet. Thanks to the preliminary information 
provided by Takagi 13 , we were able to construct a minimal 
unit-cell containing two formula units based on the hexago- 
nal structure of Na2RuOj^, a sibling compound of Na2ir03. 
The crystal structure of Na2lr03 can be viewed as an alternate 
stacking of (Ir 2 /3Na 1 / 3 )02 and Na layers. Edge-shared Ir06 
octahedra form a honeycomb lattice of Ir atoms. Na atoms 
are placed at the center of each hexagon. Upper and lower 
triangle oxygens are rotated by 3.5° to shorten the Ir-0 dis- 
tance. The positions of atoms in the unit cell were determined 
through the full structural optimization by the LDA calcula- 
tions with 0.5 xlO -3 Hatree/A of force criterion. There is a 
possible stacking disorder in the types of the Na-layers rela- 
tive to the (Ir2/3Na!/ 3 )02 layers. We have checked the ef- 
fect of different stacking sequences and observed a negligible 
change in the energy dispersions. Since the basic electronic 
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FIG. 1: (Color online) Electronic band structures of Na2lr03 within 
(a) LDA, (b) LDA+SO, and (c) LDA+U+SO schemes. Green, red, 
and blue colored energy dispersions in (a) are indicating e g , e' g , and 
aig bands respectively, induced by the largest cubic and the next 
largest trigonal crystal fields. 




structure is dominated by the in-plane Ir-0 hybridization and 
remains intact regardless of the stacking sequence, we will 
focus on the electronic structure without structural disorder 
hereafter. 

We investigated the electronic and magnetic structures of 
Na 2 Ir0 3 by performing LDA, LDA+SO, and LDA+£/+SO 
calculations. Calculated electronic band structure near E^ are 
shown in Fig.[T] The LDA band structure in Fig.[TJa) features 
the Ir 5d bands of e g and t^g components separated by a large 
cubic crystal field A cu bi C ~ 4 eV. While narrow e g bands are 
located at 3 eV above Ep, the top of t^g bands are pinned at 
Ey and spread out to -2.0 eV below Ep. Due to the extended 
nature of Ir 5d orbitals, there are large contributions to the 
band structure from both the indirect hopping via the Ir 5d-0 
2p hybridization and the direct hopping between the neigh- 
boring Ir 5d orbitals. From the tight-binding analysis 18 , even 
the next-nearest-neighbor hopping terms through oxygen and 
sodium atoms make significant contributions to the LDA band 
structure. 

The trigonal crystal field (A tr i gona i) splits the t^g bands into 
aig and e' g states. In addition, there is a strong hybridiza- 
tion between neighboring Ir 5d orbitals which gives rise to 
the bonding and anti-bonding of e' g orbitals. The bonding and 
anti-bonding doublet states consist of e' g orbital pairs of two Ir 
atoms per unit cell. At the Y point of the LDA band structure, 
the e' g anti-bonding states, to be called by cab, are close to E F 
while the e' g bonding states, to be called by e#, are at about 
—0.8 eV. The a\ Q bands located at — 1 eV have a negligible 
effect of the hybridization between neighboring Ir atoms but 
show a relatively large c-axis dispersion, which may be de- 
rived from the character of a\ Q orbitals pointing toward the 
Na atoms in the next layers. Here it is noted that the appear- 
ance of the 6ab doublet at E? in the LDA band structure of 
Na2lrC>3 is in contrast to the presence of almost degenerate 
t2 g state in Sr2lr04 which serves as a basis for the j efr =l/2 
state when SOC is introduced 9 . 

In the LDA band structure, the doubly degenerate cab 
states form a narrow band and cross Ep. The introduction 
of SOC breaks the degeneracy of cab by preserving the time- 
reversal symmetry so that the 6ab bands split off over the 
whole Brillouin zone (BZ) as shown in Fig.[TJb). Despite the 
split of eAB bands, the LDA+SO band structure is still metal- 



FIG. 2: Electronic band structures from LDA+SO calculations with 
the scaling factors of SOC strength Aso/Ao are (a) 0.5, (b) 1.0, (c) 
1.5, and (d) 2.0, where Ao is the SOI magnitude of a real Ir atom. 
Gap opens when Aso/Ao is increasing from 1.0 to 1.5. 



lie with a small electron pocket at the A point and hole pockets 
off the k c = plane near M. From the tight-binding analy- 
sis of the Na2lr03 band structure 18 , we obtained A tr i gona i ~ 
0. 6 eV , which is larger than the SOC parameter Aso ~ 0-4 
e \M12l Thus the band structure of Na2lr03 near Ep is char- 
acterized by the bonding e b and anti-bonding 6ab states with 
A cubic > A tr igonai > Aso- Since A tri gonai > Aso, however, the 
Cab character of the bands are maintained. Contrary to the 
layered perovskite Sr2lr04 system, where the SOC entangles 
almost degenerate t^g orbitals with spin states and produces 
the spin-orbit integrated j e ff=l/2, the strong trigonal field in 
Na2lr03 suppresses the mixing of a\ g and e' g states. Instead, 
the SOC acting on the e' g subspace plays a role of effective 
Zeeman coupling, the details of which will be discussed later. 
The presence of the effective Zeeman coupling is manifested 
in the parallel splitting of eAB and eg bands. 

Similarly to the case of Sr2lr04, both the on-site Coulomb 
interaction and the SOC are expected to be important in the 
description of Ir 5d states. The LDA+Z7+SO band structure 
shown in Fig.[TJc) was calculated with an effective U = 2.0 
eV, which was found to be consistent with angle-resolved pho- 
toemission and optical spectroscopy experiment^. As a result 
of the combined action of both on-site Coulomb interaction 
and SOC, a small band gap arises between the SO-split eAB 
bands. Two eAB bands form valence and conduction bands 
with nearly the same dispersion above and below Ep, respec- 
tively. Contrary to the non-magnetic metallic solution of the 
LDA and LDA+SO calculations, the LDA+/7+SO solution 
predicts an AFM ordering with local magnetic moments ly- 
ing within the ab plane. The magnitude of total moment is 
0.47 ii b per each Ir atom, which is decomposed into the spin 
moment of 0.12 /i^ and the orbital moment of 0.35 

Despite that the importance of both U and Aso, the nature 
of the insulating ground state of Na2ir03 is quite distinct from 
that of Sr2lr04. In Sr2lr04, the j e ff=l/2 degeneracy can not 
be lifted by the SOC and the Mott-Hubbard gap can be at- 
tained only when the on-site U is introduced. Thus break- 
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ing the time-reversal symmetry is essential to get the insulat- 
ing ground state of Sr2lr04. In the case of Na2lr03, how- 
ever, the broken time-reversal symmetry is not required to ac- 
quire the insulating state. As shown in Fig. [TJb), the SO- split 
6ab bands are separated over the whole BZ so that the in- 
crease of the SOC strength can enlarge the already present 
gap between two cab bands. To probe this idea, we carried 
out DFT calculations by controlling the SOC strength, which 
can be achieved by changing the scaling factor when generat- 
ing the j-dependent pseudo-potential 14 . Calculated results for 
the scale factors Aso/Ao= 0.5, 1.0, 1.5, and 2.0 are shown in 
Fig. [2] Taking the original SOC in the real Ir atom as Ao as 
a reference, Aso/Ao= 1.5 was found to be enough to open a 
full insulating gap. We call these insulating ground states as 
spin-orbit (SO) insulators, which have energy gaps generated 
by the SOC. SO insulators have no local moment and pre- 
serve the time-reversal symmetry and thus are distinct from 
the Mott-Hubbard insulator. 

To understand the origin of SO insulators, we consider the 
SOC matrix elements within the e' g subspace. Since the de- 
generate e' g states can be written by 



□ : Mott insulator 
a : Ferromagnetic metal 
o : Paramagnetic metal 
• : Spin-orbit insulator 
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FIG. 3: Phase diagram in the \so-U parameter space depicting four 
different phases from LDA+[/+SO calculations with varying U and 
Aso values. Paramagnetic metallic phase appears near the origin, 
Mott insulator in the region of U > 1, and SO insulator in the region 
of Aso > 1 and U < 1. The real ground state is located inside Mott 
insulating territory. 
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where the basis sets are | e' )<8> \ S = \) = {| e[a), \ e[/3), \ 
e f 2 a), | e 2 (3)} and h is the unit vector along the c-axis, i.e., 
the [1 1 1] direction in the local coordinate of Ir06 octahedron. 
This block-diagonal form comes from the fact that (L) is si- 
multaneously diagonalized within e' g manifold and its eigen- 
values are n and — h, respectively. Here the SOC terms in e' g 
did as an internal magnetic field perpendicular to the afr-plane. 
The internal field gives rise to an effective Zeeman splitting, 
but the field direction in the e[ component is opposite to that 
in the e' 2 component. Thus, the effective Zeeman coupling 
does not break the time-reversal symmetry and \e[a)-\e 2 P) 
and \e f 1 /3}-\e f 2 a) remain as time-reversal partners. Since the 
eAB states are the anti-bonding combination of the e' g orbitals 
of neighboring Ir atoms, the splitting of eAB bands by the ef- 
fective Zeeman coupling is proportional to the SOC strength 
as shown in Fig.[2| especially at the T point. 

The ground states of Na2ir03 with the large SOC strength 
are SO insulators. The band gaps are induced by the effective 
Zeeman coupling of the SOC within the e' g subspace. Their 
characters are different from other types of band insulators 
such as covalent or ionic ones. The Fermi level is placed be- 
tween bonding and anti-bonding bands in covalent solids and 
between different ionic configurations in ionic solids. In SO 
insulators, the gap is not driven by bonding characters, but 
mainly related to the symmetry of the states at E F . In a sense 
that their band gaps are generated by the SOC, SO insulators 



share the same ground with topological insulators though it is 
necessary to prove the non-trivial topology of its ground state. 

One important consequence of the SO insulating phase is 
the proximity of the AFM ground state to the SO insulator 
state. In the LDA+Z7+SO calculation, the AFM ordered lo- 
cal moments are aligned in the afr-plane. Due to the huge 
internal field along [111] direction, it is hard to break the 
time-reversal symmetry and to develop local magnetic mo- 
ments along that direction. Thus, transverse magnetic mo- 
ments which are perpendicular to the internal field can be eas- 
ily developed. Strong magnetic anisotropy originated from 
the internal magnetic fields might be seen in magnetic suscep- 
tibility measurements. 

To elucidate the relation between SO insulator and AFM 
Mott insulator phases, we explored a possible phase dia- 
gram of Na2lr03 in an extended \so~U parameter space and 
present the result in Fig. [3] When U is small and Aso / Ao 
is less than 1.5, the ground state remains as a paramagnetic 
metal. When there is no SOC, i.e., Aso = 0, a ferromagen- 
tic metallic phase develops in a narrow range of the parameter 
space with Aso = upto U = 5.0 eV. This ferromagnetic 
state becomes unstable in the presence of the SOC. On the 
other hand, for the value of U smaller than about 1 eV, the SO 
insulator phase emerges as a non-magnetic insulator. Since 
the band gap is induced by the effective Zeeman coupling of 
the SOC within the e' g subspace, the Kramers degeneracy of 
the valence states holds up and the time-reversal symmetry 
remains unbroken. For the finite Aso> Mott insulating AFM 
states develop as U becomes larger than about 1.0 eV. The 
difference between two insulating phases, i.e., the criterion 
for the boundary is the existence of local magnetic moments. 
The Mott insulating phase has AFM ordering where on-site 
Coulomb repulsion breaks the symmetry developing local mo- 
ments during the correlation gap opens. Our LDA+Z7+SO 
calculation predicts that the real ground state of Na2ir03 is 
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FIG. 4: (Color online) Schematic drawing of the (Ir 2 /3Nai/ 3 )02 
plane and magnetic configuration of the AFM insulating ground state 
of Na2lr03- Magnetic moments are ordered anti-ferromagnetically 
lying on the afc-plane due to the strong internal field along the c-axis. 
Not only the NN exchange J (dotted arrow) but the NNN exchange 
J' (dashed arrow) are significant and may give rise to magnetic frus- 
tration. 



where G is the one-particle Green's function and V is on-site 
exchange interaction potential 20 . The result is J' / J = 0.47, 
which means that NN and NNN exchange coupling strength 
are comparable and they might be a source of frustration. 
Above result is mainly attributed to the extended nature of Ir 
5d orbitals. Large direct overlap between NN Ir atoms gives 
FM direct exchange interaction, competing with AFM su- 
perexchange from oxygen mediated hopping channels and fi- 
nally reducing AFM exchange J. On the other hand, the NNN 
hopping integrals are not negligible that the NNN AFM inter- 
action J' can be comparable and frustrate long range AFM 
ordering. 

In conclusion we have shown that the spin-orbit entangled 
e' g states under the strong internal Zeeman field driven by the 
SOC lead to an unusual band gap. The predicted AFM ground 
state is in close proximity to the SO insulator phase where the 
AFM ordering in Na2lr03 becomes strongly anisotropic with 
quenched moments along the c-axis. The highly anisotropic 
AFM state in Na2lr03 may serve as a model system for the 
two-dimensional XY model with frustrated exchange interac- 
tions. One may be able to drive a crossover between AFM and 
SO insulators through the modulation of structural parameters 
or chemical substitution, though we need more study on the 
role of SOC in the Mott AFM phase in connection with the 
topological nature of SO insulators. 



located in the Mott AFM region with U = 2.0 ~ 3.0 eV and 
AsoMo = 1- However, the strongly anisotropic nature of its 
AFM ordering originates from its proximity to the SO insula- 
tor phase. 

Another important aspect in Na2lr03 is magnetic frus- 
tration indicated in large #cw/^n ratio from susceptibility 
measurements-^. To reveal the origin of frustration, we have 
estimated exchange interactions J and J' between nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) Ir atoms re- 
spectively. (Fig. [4]) Calculation scheme is based on the pertur- 
bation formalism. 
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